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Correlation  From  Randomness:  Quantitative  Analysis  of  Sputtered  Graphite 
Surfaces  Using  the  Scanning  Tunneling  Microscope 

Elliott  A.  Eklund*l*3^  Eric  j  Snyder2i3  and  R.  Stanley  Williams^.S 

^Department  of  Physics 
^Department  of  Chemistry  and  Biochemistry 
SSolid  State  Science  Center 
University  of  California  Los  Angeles,  CA  90024 

Abstract:  The  Scanning  Tunneling  Microscope  (STM)  was  used  to  quantitatively 
examine  the  surface  morpholog>'  of  highly  oriented  pyrolytic  graphite  (HOPG) 
surfaces  bombarded  with  5  keV  Ar'*’  ions.  Constant  current  topographs  clearly 
showed  that  the  morphology  of  the  resulting  nonequilibrium  surfaces  depended 
sensitively  on  the  ion  flux,  the  ion  fluence  and  the  sample  temperature.  For  low 
and  intermediate  ion  fluence,  an  increase  in  surface  roughness  with  increasing 
ion  flux  w’as  observed.  For  constant  ion  flux,  the  surfaces  developed  structures 
with  heights  proportional  to  the  ion  fluence  and  separations  characterized  by  a 
correlation  length  that  diverged  with  increasing  fluence.  Increasing  the  sample 
temperature  during  bombardment  produced  smoother  surfaces  as  a  result  of 
enhanced  surface  diffusion.  The  autocovariance  function  G(L)  and  the  height 
correlation  function  <lh(q)|2)in  reciprocal  space  were  calculated  directly  from 
the  STM  topographs.  The  latter  was  then  compared  with  a  linear  response  theoiy 
for  the  formation  of  rough  surfaces  under  nonequilibrium  conditions  and  a 
scaling  analysis.  The  fluence  dependence  of  the  surface  roughening  of  graphite 
could  not  be  explained  by  shot  noise  alone  in  a  linear  theory,  whereas  surface 
diffusion  and  redeposition  of  sputtered  material  satisfactorily  accounted  for  the 
correlation  observed  in  the  bombardment-induced  features. 


I.  Introduction 


One  of  the  fundamental  challenges  in  materials  science  is  to  understand 
the  effects  of  particle  radiation  on  the  nature  of  solid  surfaces. For  example, 
the  properties  of  sputter  deposited  thin  films  depend  sensitively  on  the  surface 
topology,  which  is  determined  by  the  sputter-growth  process.^  A  similar 
situation  exists  for  ion-beam  sputter  etching,  a  widely  used  technique  in  surface 
science  and  microelectronics,  which  also  produces  nonequilibrium  surfaces  with 
topography  that  depends  on  the  sputtering  conditions. ^ 

Erosion  of  surfaces  by  sputtering  is  important  in  a  great  variety  of 
disciplines.  In  surface  science  experiments,  clean  surfaces  are  usually  prepared 
by  repeated  cycles  of  sputtering  with  noble  gas  ions  and  annealing  at  elevated 
temperatures.®  In  depth-profiling  analyses,  sputter  etching  is  used  to  strip  away 
successive  layers  of  the  sample,  exposing  sub-surface  material  for  compositional 
determination. 2  With  this  technique,  control  of  the  surface  roughness  is 
important  because  the  bombardment  induced  topography  determines  the  depth 
resolution  of  the  compositional  analysis.'^*®  Focused  ion  beam  raicromachining 
has  been  used  in  the  electronics  industry  for  the  repair  of  X-ray  lithography 
masks  and  the  fabrication  of  device  microstructures.®  In  some  cases,  sputtering 
is  an  imwanted  but  important  side  effect  of  the  environment,  as  in  the  erosion  of 
the  inner  walls  of  fusion  plasma  containment  vessels  or  the  outer  skins  of  space 
vehicles  as  they  pass  through  the  atmosphere.  In  fact,  fusion  plasma  chambers 
commonly  have  a  graphite  interior  lining^®  and  carbon  composites  are  highly 
desirable  materials  for  space  vehicles  because  of  their  strength  and  light 
weight. The  experiments  presented  in  this  study  are  particularly  relevant  to 
the  understanding  of  the  sputtering  of  graphite  and  carbon  composites,  since  ion  -- 
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bombardment  shoiild  produce  similar  amorphous  surfaces  on  both  types 
materials  after  a  short  time. 

In  general,  surfaces  exposed  to  high  levels  of  particle  radiation  develop 
characteristic  morphologies.  On  the  macroscopic  level  (I'.e.,  length  scales  larger 
than  one  micron),  the  effects  of  particle  radiation  are  well  understood.  Ion 
bombardment  erodes  edges,  comers  and  peaks  preferentially  because  these 
prominent  features  are  exposed  to  a  larger  incident  flux.  This  radiation  erosion 
is  further  enhanced  by  surface  diffusion  and  evaporation/recondensation.  In  the 
continuum  description  of  this  process,  one  either  assumes  that  the  local  erosion 
rates  are  proportional  to  the  local  incoming  flux,  or  one  uses  the  so-called 
"Huygens  Construction", in  which  successive  profiles  of  the  eroding  surface 
are  treated  as  a  propagating  kinetic  wave.  In  this  latter  method,  one  calculates 
the  surface  profile  produced  after  a  time  At  by  constructing  spherical  erosion 
zones,  which  have  radii  proportional  to  the  ion  flux,  the  local  sputtering  yield  and 
the  time  step  At,  at  each  point  on  the  starting  surface  profile.  The  new  surface 
profile  is  then  given  at  each  point  by  the  envelope  of  these  spherical  erosion 
volumes.  Numerical  simulations  of  the  macroscopic  evolution  of  sputtered 
surfaces  have  produced  results  that  compare  favorably  with  experiments. 

1  fi 

At  the  sub-micron  level,  our  understanding  is  much  less  complete. 
Electron-microscope  studies  of  sputtered  surfaces  indicate  that  particle  radiation 
will  roughen  the  surface  rather  than  smooth  it,  and  prominent  cone-like  features 
have  been  reported.  The  stochasticity  of  the  bombardment  process  can  be 
included  in  theoretical  models  by  adding  a  term  describing  shot  noise,  the 
variation  of  the  incident  flux  around  some  mean  value.  Both  analytical  studies 
and  Monte  Carlo  simulations20.21  of  the  effect  of  shot  noise  on  continuum 
theories  of  growth  by  sputter  deposition  indicate  that  such  nonequilibrium 
surfaces  evolve  roughness.  The  features  produced  are  characterized  by  a 
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correlation  length  ^(t)  which  diverges  with  deposition  time  t.  Similar  results  are 
expected  for  the  case  of  sputter  etching,  but  until  now  there  has  been  little 
experimental  confirmation  of  this  description. 

Our  use  of  the  Scanning  Txinneling  Microscope  (STM)22,23  to  examine  the 
surface  topography  of  graphite  after  sputter  etching  offered  many  unique 
opportunities  for  the  study  of  nonequilibrium  surfaces.  Unlike  the  scanning 
electron  microscope  (SEM),  which  has  been  the  primary  analysis  tool  for  imaging 
the  topography  of  sputtered  materials,  the  STM  can  quantitatively  measure  the 
height  profile  h(r)  of  the  surface  over  a  l^'ide  range  of  lateral  and  vertical  length 
scales,  with  a  lateral  resolution  of  better  than  10  A  and  a  vertical  resolution  of  a 
fraction  of  an  Angstrom.22,23  This  allowed  us  to  calculate  from  the  topographs 
the  auto-covariance  function, 24  defined  as; 

G(lr2-rjl)  =  (h(r2)h(rj)>  -  (h(r))2  ,  (1) 

where  h(r- )  is  the  height  of  the  surface  at  position  r-  and  L  =  I  r2  -  I  is  a 
separation  length  on  the  surface,  and  the  height  correlation  function,  defined  as: 

<lh(q)|2)  =  eiq-r  <(h(0)-h(r))2)^/Area,  (2) 

with  (•••  indicating  a  sample  average  after  t  seconds  of  exposure.  The 
calculation  of  the  auto-covariance  function  G  provided  a  quantitative 
determination  in  real  space  of  the  base  to  peak  height  variation  and  the  short 
range  lateral  correlation  of  the  surface  features.  Calculation  of  the  height 
correlation  function  (lh(q)  1 2)  allowed  us  to  determine  the  reciprocal  space 
scaling  behavior  of  the  surfaces  as  a  function  of  the  wavenumber  q  and  to  make 


4 


comparisons  of  our  experimental  data  to  existing  theories  of  surface  morphology 
created  by  sputter  erosion.  In  other  surface  roughening  problems,  e.g. ,  thermal 
roughening,  this  correlation  function  was  found  to  be  a  good  measure  of  the 
surface  geometry  .25  Because  of  loss  of  phase  information  in  the  height 
correlation  function,  however,  it  is  less  sensitive  to  uncorrelated  -  but  prominent  - 
-  siirface  features,  which  are  much  more  easily  visible  in  the  real  space  STM 
topographs. 

The  STM  has  been  used  by  other  researchers  to  study  rough  surfaces 
produced  by  fracture  and  deposition  by  sputtering  and  evaporation.26'28  These 
studies  have  addressed  some  of  the  special  problems  associated  ^^ith  the  scanning 
of  rough  surfaces  with  a  tunneling  microscope,  in  particular,  that  the  STM 
topographs  can  depend  on  the  shape  and  sharpness  of  the  tip  and  on  the 
roughness  of  the  surfaces.  Asimimetric  tips  will  produce  artificial  anisotropy  in 
the  images  while  blunt  tips  will  fail  to  reproduce  small,  closely  spaced  features.29 
Surfaces  with  very  sharp  protrusions  will  yield  images  of  the  sides  of  the 
tunneling  tip  instead  of  the  surface  features.30  We  have  been  particularly  careful 
to  check  that  tip  related  artifacts  have  not  contributed  seriously  to  the  topographs 
in  this  study.  Other  factors  that  must  be  considered  include  the  maximum  size  of 
the  tunneling  microscope  scans  and  the  density  of  the  data  points  in  a  topograph. 

A  major  difference  between  this  and  prerious  studies  of  sputter-induced 
topography  was  that  the  freshly  cleaved  graphite  surfaces  were  atomically  flat 
over  tens  of  thousands  of  square  Angstroms^^  before  bombardment.  This  allowed 
us  to  investigate  the  initial  stages  of  sputter  induced  morphology  from  flat 
surfaces  starting  with  low’  ion  fluences  (10l6ions/cm2)  and  to  characterize  the 
evolution  of  the  surface  profile  as  a  function  of  time.  The  real  space  topographs, 
the  auto-covariance  functions  and  the  height  correlation  functions  provided 
complementary  information  that  allowed  us  to  investigate  the  dependence  of  the 
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s\irface  morphology  of  sputtered  graphite  on  several  of  the  most  commonly  varied 
sputtering  parameters.  The  surfaces  were  sensitive  to  variations  in  ion  flux,  ion 
fluence  and  sample  temperature.  In  particular,  for  the  lowest  ion  fluences,  the 
base  to  peak  height  variations  of  the  sputtered  surfaces  increased  with  increasing 
ion  flux.  For  all  fluxes  used,  the  correlation  length  and  the  vertical  size  of  the 
surface  features  increased  as  the  fluence  was  increased.  Moreover,  the  size  of  the 
features  scaled  approximately  linearly  with  the  fluence,  which  is  surprising  for 
an  essentially  stochastic  process.  Sputtering  at  elevated  sample  temperatures 
produced  surfaces  that  were  smoothed  by  enhanced  surface  diffusion. 

II.  Experimental  Procedure 

The  STM  used  in  these  experiments  was  of  the  flip  stage  design  introduced 
by  Demuth  et  al.32  Our  instrument  was  designed  to  operate  at  atmospheric 
pressure  and  a  variety  of  precautions  were  taken  to  make  the  microscope 
insensitive  to  mechanical,  electrical,  and  thermal  disturbances  in  the  laboratory. 
The  mechanical  isolation  of  the  instrument  w’as  accomplished  in  two  stages;  the 
STM  assembly  was  placed  on  a  platform  hung  from  rubber  cords  to  isolate  it  from 
low  frequency  room  vibrations  (1-10  Hz)  and  the  microscope  stage  itself  rested  on 
another  isolation  stage  of  steel  plates  separated  by  viton  rubber^S  which  served  to 
damp  out  higher  frequencies.  Isolation  of  the  STM  electronics  was  accomplished 
by  careful  shielding  of  sensitive  components  and  by  making  all  ground 
cormections  to  a  single,  high  quality  ground  point.34  For  the  experiments  in  this 
study,  the  total  electronic  and  mechanical  noise  was  less  than  2%  of  the  tunneling 
signal,  which  meant  that  the  rough  surface  features  seen  in  the  STM  topographs 
of  sputtered  graphite  had  no  significant  system  noise  contributions.  Thermal 
drift  is  a  serious  problem  in  many  tunneling  microscopes  operating  at 
atmospheric  pressure  and  can  cause  distortion  in  STM  images.  The  constant 
current  imaging  mode  w'ith  its  slow^er  scanning  speeds  was  needed  for  the 
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roughened  graphite  samples  and  required  approximately  2  minutes  to  acquire  a 
full  topograph,  making  the  images  susceptible  to  drift.  The  thermally  induced 
drift  rates  for  our  instrument  were  low,  however,  approximately  1  A/min  for  the 
X  and  Y  scanning  axes,  and  had  a  negligible  effect  on  the  images.  A  final  level  of 
acoustical,  electrical  and  thermal  isolation  was  provided  by  an  enclosure  lined 
with  acoustically  absorbant  foam  and  a  grounded  copper  mesh  which  surrounded 
the  STM  during  scanning. 

The  3  axis  piezo-electric  drive22  that  positioned  the  tunneling  tip  was 
carefully  calibrated  before  any  distance  measurements  were  made  on  samples 
with  disordered  topography.  Two  independent  methods  were  used  to  calibrate  the 
STM.  The  atomically  resolved  features^S  seen  on  graphite  at  high  magnification 
were  used  to  calibrate  the  piezo  drives  for  small  displacements  (image  size  of  25  A 
X  25  A).  This  calibration  was  supplemented  v^nth  interferometric 
measurements^®  that  required  the  piezo-electric  elements  to  be  scanned  several 
thousand  Angstroms  at  a  time,  resulting  in  sensitirities  for  the  X,  Y  and  Z  axes  of 
12.91  ±  0.33,  12.28  ±  0.34  and  9.25  ±  0.29  A/Volt,  respectively.  The  tunneling  tips 
were  easily  prepared  by  cutting  PlTlh  thermocouple  wire,  chosen  because  it  does 
not  form  an  insulating  oxide  layer  in  air,  at  an  acute  angle  with  ordinary  wire 
cutters.  It  has  been  shown  that  mechanical  cutting  produces  tips  with  many- 
jagged  protrusions,  one  of  w'hich  acts  as  the  actual  tunneling  tip. 37  This 
technique  of  tip  preparation  has  been  used  successfully  by  other  groups38.39  and 
routinely  yielded  atomic  resolution  in  our  imaging  of  clean,  unbombarded 
graphite. 

The  surface  chosen  for  this  investigation  was  the  (0001)  face  of  highly 
oriented  pyTolytic  graphite  (HOPG),  w-hich  was  easily  cleaved  with  adhesive  tape 
to  produce  large,  atomically  flat  regions. 26  This  surface  was  inert  in  air  and 
easily  imaged  with  the  STM.3®  Graphite  also  has  a  rigid  lattice,  as  seen  by  its 
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melting  temperature  of  -3800®C.  This  indicates  that  surface  diffusion  effects  were 
minima]  and  that  topographical  features  produced  by  the  bombardment  were 
'frozen  in’  and  could  be  observed  wdth  the  STM  long  after  the  sputtering  had  taken 
place.  In  general,  surface  diffusion  is  an  important  effect  which  smooths  out 
surface  features,  so  materials  that  can  anneal  at  ambient  temperatures  are 
unsuitable  for  the  experiments  described  here. 

The  constant  current  or  topographic  mode22  was  used  to  collect  the  STM 
images,  with  a  tunneling  current  of  0.5  nA  and  a  sample  to  tip  bias  of  -100  mV. 
Topographs  collected  at  the  largest  possible  scan  size  (2400  A  x  2400  A)  show'ed 
large,  atomically  flat  areas,  while  at  higher  magnifications  (eg.  25  A  x  25  A)  the 
familiar  atomic  scale  features  of  clean  graphite  were  easily  observed. 35 

After  stable  images  of  clean  graphite  were  obtained  at  both  high  and  low 
magnifications,  the  samples  were  transferred  to  the  sample  treatment  chamber 
of  a  KRATOS  XS.4M-800  surface  analytical  system.  The  graphite  surfaces  w'ere 
then  sputter  etched  at  room  temperature  v,ith  a  beam  of  5  keV  Ar**"  ions,  rastered 
over  a  9  mm^  area  on  the  sample  and  incident  at  an  angle  of  60'  to  the  surface 
normal.  The  flux  incident  on  the  sample  was  determined  by  using  an 
electrometer  to  measure  the  ion  beam  current.  A  small  positive  bias  (45  volts)  was 
applied  to  the  sample  to  suppress  secondar>’  electron  emission.  The  experimental 
parameters  that  were  varied  in  this  study  were  the  flux  J,  the  ion  fluence  Q^Jt, 
and  the  sample  temperature  T. 

This  paper  highlights  the  results  obtained  from  over  1,000  STM  topographs 
collected  from  more  than  15  bombarded  graphite  samples.  The  three  ion  fluxes 
reported  here  were  Ji  =  6.9  x  10^3  ^  J2  =  3.5  x  lO^^  and  J3  =  6.9  x  10^^  ions/cm2-sec, 
corresponding  to  ion  beam  currents  of  approximately  1,  5  and  10  4A,  respectively. 
The  flux  was  varied  by  using  different  flow  rates  of  Ar  gas  in  the  ion  gun.  By 
varying  the  time  of  exposure  to  the  ion  beam,  the  total  fluences  obtained  w’ere  Qi  = 
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10^6,  Q2  =  and  Q3  =  10^8  ions/cm2.  A  sputtering  yield  of  1.5  carbon  atoms 
ejected  per  incident  argon  ion  was  estimated  by  Monte  Carlo  simxilations  (TRIM 
code)^0  and,  with  3.8  x  10^8  carbon  atoms/cm2  the  graphite  surface,  the  ion 
fluences  given  above  corresponded  to  the  removal  of  approximately  4,  40  and  400 
monolayers  of  material,  respectively,  from  the  graphite  surfaces.  To  investigate 
the  effect  of  temperature  (ie.  enhanced  surface  diffusion)  on  sputter  induced 
morphology,  etching  was  also  performed  on  samples  heated  to  approximately 
600K  and  900K  using  the  intermediate  flux  and  fluence  (J2  and  Q2  ) .  The  samples 
used  for  the  variable  temperature  experiments  were  re-cleaved  before  each 
bombardment  to  return  to  the  initial  flat  surface  conditions.  In  addition,  one 
graphite  sample  was  annealed  at  600  K  after  clearing  to  release  any  surface  stress 
and  was  allowed  to  cool  before  bombardment.  The  surfaces  produced  were 
similar  to  those  sputtered  without  annealing  before  sputtering,  indicating  that  the 
effects  of  clearing  on  surface  morpholog>’  created  by  ion  bombardment  were 
minimal. 

The  graphite  samples  were  re-examined  with  the  STM  after  etching  using 
identical  operating  parameters  and.  if  possible,  the  same  tunneling  tip  used  prior 
to  sputtering.  The  results  were  reproducible  from  sample  to  sample  and  even 
with  different  tips,  indicating  that  tip  artifacts  did  not  seriously  affect  the  data. 
Topographs  with  serious  artifacts  were  observed  but  were  discarded  from  the  data 
set.  Fourteen  different  image  sizes,  ranging  from  2400  A  x  2400  A  to  25  A  x  25  A, 
were  used  to  investigate  each  of  the  sputtered  surfaces.  Typically  over  50  images  of 
various  sizes  were  collected  at  each  location  and  several  different  locations  were 
imaged  on  each  sample.  The  topographical  data  for  each  image  were  digitally 
stored  by  the  STM  computer  as  arrays  of  118  x  118  data  points  for  subsequent 
riewing  and  calculations.  In  addition,  the  images  were  leveled  by  a  least  squares 
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fit  of  the  data  to  a  plane,  thereby  removing  any  effects  of  sample  tilt  but  also 
possibly  remo\'ing  long  wavelength  features  from  the  data. 

The  bombardment  by  5  keV  argon  ions  was  expected  to  disrupt  the  graphite 
surfaces  and  possibly  make  them  more  prone  than  freshly  cleaved  samples  to 
reaction  with  air.  In  fact,  after  2  -  3  days  of  scanning  sputtered  surfaces  with  the 
STM,  the  images  became  noisy  and  some  of  the  reproducibility  was  lost.  These 
problems  were  avoided  for  the  samples  used  in  this  study  by  minimizing  the 
exposure  of  the  surfaces  to  air  during  the  time  needed  to  acquire  the  topographs. 
At  all  other  times,  the  samples  were  kept  under  vacuum  in  the  etching  chamber 
or  were  stored  in  a  nitrogen  filled  glovebox,  limiting  the  total  air  exposure  for  each 
surface  to  a  maximum  of  a  few  hours. 

III.  Results  and  Qualitative  Discussion  of  Topographs 

The  topographs  presented  in  Figs.  1  -  4  show  the  beha\nor  of  the  graphite 

c  • 

surfaces  at  an  image  size  of  2400  A  x  2400A  as  the  flux  J,  fluence  Q  and 
temperature  T  are  varied.  Each  of  the  Figs.  1-3  shows  a  different  ion  fluence, 
and  the  effect  of  changing  the  ion  flux  for  fixed  fluence  is  shown  within  each  of 
these  figures.  Fig.  4  illustrates  the  effect  of  sputtering  at  different  temperatures. 
The  reproducibility  and  scaling  of  surface  features  as  the  magnification  is 
increased  is  demonstrated  in  Figs.  5  -  7. 

The  dependence  of  the  morphology  on  the  etching  flux,  for  fixed  fluence,  is 
clearly  seen  in  Figs.  1  -  3.  For  the  lower  fluences,  and  Q2,  Figs.  1  and  2  show 
that  the  surfaces  became  steadily  rougher  as  the  ion  flux  was  increased.  The 
three  fluxes  J^,  J2  and  J3  were  in  the  ratio  of  ~1:5;10.  The  features  in  Figs,  la,  lb 
and  Ic  have  heights  of  -4 , 20  and  40  A,  respectively,  while  those  in  Figs.  2a,  2b 
and  2c  have  heights  of  ~10,  30  and  60  A,  roughly  proportional  to  the  ratio  of  the 
fluxes.  Also  apparent  in  Figs.  1  and  2  is  the  increasing  amplitude  of  the  short 
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wavelength  roughness  as  the  flux  was  increased.  For  the  lowest  flux,  Figa.  la 
and  2a  show  features  separated  by  -200  A  while  Figs.  Ic  and  2c,  showing  surfaces 
sputtered  with  the  highest  flux,  are  dominated  by  features  with  separations  <100 
A.  Featxires  with  both  long  and  short  length  scales  are  seen  for  surfaces 
sputtered  at  the  intermediate  flux,  shown  in  Figs,  lb  and  2b.  These  trends  are 
modified  for  the  highest  fluence  Q^,  shown  in  Fig.  3.  The  separation  of  surface 
features  did  not  change  significantly  between  the  two  lower  fluxes,  since  Figs.  3a 
and  3b  both  show  similar  large  structures  separated  by  -1500  -  2000  A,  but  did 
decrease  for  the  sample  sputtered  at  the  highest  flux,  which  shows  prominent 
features  at  smaller  separations  of  -800  A  in  Fig.  3c. 

The  most  striking  effect  seen  in  comparing  Figs.  1  -  3  is  the  sensithity  of  the 
surface  morphology  to  the  fluence  Q,  for  all  fluxes.  In  general,  at  the  fluence  of 
lO^®  ions'cm2  shown  in  Fig.  1,  the  surfaces  all  exhibited  short  range  roughness 
of  small  amplitude,  with  features  4  -  40  A  in  height  separated  by  distances  of  50  ■ 
200  A.  Increasing  Q  to  10 ions/cm^  produced  surfaces  dominated  by  much 
larger  features,  60  -  200  A  high,  with  separation  distances  that  approached  the 
image  size,  as  seen  in  Fig.  3.  The  transition  between  these  extremes  is  shown 
clearly  in  Fig.  2  for  surfaces  sputtered  at  the  intermediate  fluence  of  10^"^ 
ions'cm^  ,  where  4  -  10  A  high  features  at  small  separations  are  seen  along  with 
hillocks  that  are  2  to  3  times  higher  at  separations  of  several  hundred  Angstroms. 

Because  of  enhanced  surface  diffusion,  it  is  expected  that  sputter  etching  of 
samples  at  elevated  temperatures  should  yield  surfaces  that  are  less  rough  than 
those  which  are  sputtered  at  ambient  conditions.  This  prediction  is  confirmed  in 
Figs.  4a,  4b  and  4c,  which  show  graphite  surfaces  bombarded  at  temperatures  of 
300K  (ambient),  600K  and  900K,  respectively,  using  the  intermediate  flux  J2  and 
fluence  Q2.  From  Fig.  4b  we  see  that  raising  the  sample  temperature  T  to  600K 
produced  surfaces  which  had  roughly  the  same  heights  (10  -  20  A)  as  those 
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produced  at  ambient  conditions,  but  the  characteristic  wavelengths  of  the  features 
were  longer.  Fig.  4c  shows  that  a  further  increase  in  sample  temperature  to  T  = 
900K  resulted  in  significantly  smoother  surfaces.  For  these  conditions,  the 
heights  of  the  features  were  observed  to  be  only  3  •  5  A. 

As  mentioned  above,  topographs  were  acquired  at  many  different  image 
sizes  for  each  surface.  This  is  important  for  several  reasons.  To  eliminate  the 
possibilities  of  tip  artifacts  and  problems  with  mechanical  noise,  it  is  critical  to 
make  certain  that  the  features  observed  in  larger  images  scale  properly  as  the 
scan  size  is  reduced.  In  addition,  the  STM  is  able  to  probe  the  small  scale 
structures  more  acciorately  at  the  smaller  scan  sizes.  This  is  a  result  of  the 
maximum  data  density  of  118  points  per  scan  line  and  of  the  fact  that,  for  constant 
scanning  frequency,  the  tip  velocity  decreases  as  the  image  size  is  reduced, 
allo\sing  more  time  for  the  STM  to  respond  to  the  the  surface  topography.  WTiile  it 
is  impractical  to  present  all  topographs  taken  at  different  sizes,  it  is  useful  to 
examine  a  representative  sampling.  Therefore,  images  taken  at  scan  sizes  of  1000 
A  X  1000  A,  250  A  X  250  A  and  50  A  x  50  A  are  shown  in  Figs.  5,  6  and  7  for  the 
surfaces  imaged  at  a  size  of  2400  A  x  2400  A  in  Figs,  la,  2c  and  3b,  respectively. 
These  surfaces  represent  the  limits  of  roughness  in  the  data  set:  one  relatively 
smooth  (Fig.  la),  one  with  large  amplitude,  short  wavelength  features  (Fig.  2c) 
and  one  with  large  structures  at  large  separations  (Fig.  3b).  The  structures 
imaged  at  low  magnification  (2400  A  x  2400  A)  in  Figs.  1  -  3  are  seen  to  scale 
properly  dowm  to  the  250  A  x  250  A  image  size  in  figs.  5*7,  confirming  the  size 
measurements  of  these  features.  In  addition,  if  the  scale  for  the  image  was 
reduced  enough,  the  surfaces  eventually  appeared  to  be  smooth  (note  that  all 
topographs  have  been  leveled,  so  a  region  that  was  on  the  side  of  a  hillock  is 
showm  as  flat),  as  observed  in  Figs.  5c,  6c  and  7c  at  an  image  size  of  50  A  x  50  A. 
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In  summary,  it  is  apparent  from  this  set  of  STM  images  that  the 
noneqmlibrium  surfaces  resulting  from  the  sputtering  of  cleaved  graphite 
substrates  are  strongly  dependent  on  the  bombardment  conditions.  For  the  lower 
fluences,  the  height  of  the  features  on  the  sputtered  surfaces  increased 
approximately  linearly  with  increasing  flux  J.  In  addition,  for  each  ion  flux 
investigated,  the  surface  morphology  depended  strongly  on  fluence  Q.  In  general, 
short  wavelength  features  were  produced  at  the  lowest  fluence  while  the  features 
increased  in  amplitude  and  separation  as  Q  was  increased.  Sputter  etching 
produced  progressively  smoother  surfaces  with  increasing  sample  temperature. 
Finally,  an  investigation  of  images  of  bombarded  graphite  taken  at  various 
magnifications  indicated  that  surface  features  scaled  correctly  as  the  image  size 
was  reduced,  confirming  the  reliability  of  the  data  at  large  and  small  scan  sizes. 
This  scaling  investigation  also  revealed  that  all  of  the  sputtered  surfaces  began  to 
appear  smooth  at  an  image  size  of  -  50  x  50  A^. 

rV.  Quantitative  Analyses 

The  STM  topographs  clearly  show  trends  in  surface  morpholog>'  as  the 
sputtering  parameters  are  varied  and  allow  qualitative  comparisons  of  surface 
roughness  between  the  samples.  The  value  of  tunneling  microscopy  for  this 
study,  however,  lies  in  the  ability  to  calculate  from  the  height  data  in  the 
topographs  the  statistical  functions  that  quantitatively  describe  the  scaling 
behavior  of  the  surface  profiles.  The  quantities  calculated  for  the  graphite 
surfaces  in  this  study  are  the  corrugation,  the  auto-covariance  function  and  the 
height  correlation  function. 

The  corrugation  o  is  the  standard  deviation  of  the  surface  slope  and  is 
calculated  as  a  function  of  lateral  length  for  distances  of  2  -  1009r  of  the  image 
size.  The  corrugation  values  are  used  to  check  for  tip  imaging  artifacts  in  the 
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topographs.  A  large  value  of  o  at  a  partictilar  length  would  indicate  that,  at  that 
separation,  the  protrusions  and  depressions  on  the  surface  have  aspect  ratios 
similar  to  the  tunneling  tip  and  that  they  can  effectively  act  to  image  topography 
on  the  sides  of  the  tip.  This  would  produce  an  image  which  is  a  complicated 
convolution  of  the  tip  shape  and  the  actual  surface  profile.  Small  values  of  the 
corrugation  confirm  that  the  tip  is  shaip  compared  to  surface  details  and  that  the 
STM  is  able  to  probe  the  features  correctly,  producing  an  image  that  accurately 
represents  the  true  topography.  As  an  example,  consider  the  saw*tooth  shaped 
surface  profile  with  features  of  height  L  shown  in  Fig.  8.  At  the  lateral  length 
scale  of  L,  the  slopes  of  the  sides  of  the  features  are  +1  and  -1,  the  average  slope  is 
zero  and  the  corrugation  o  is  1.  The  apex  angle  of  each  feature  is  90°,  indicating 
that  for  o  =  1,  a  reasonable  tunneling  tip  should  be  able  to  accurately  probe  this 
artificial  surface  profile  without  serious  convolution  of  the  tip  shape  into  the 
image,  except  at  the  concave  comers.  Our  experience  has  shown  that  o  values 
less  than  1.0  do  not  produce  serious  tip  artifacts. 

The  topographs  of  rough  or  random  surfaces  contain  an  overabundance  of 
information,  making  it  difficult  to  pick  out  important  trends  and  subtle  features. 
To  aid  in  the  data  analysis,  we  have  calculated  from  the  STM  topographs  the  auto¬ 
covariance  functions  G(L),  given  by  Eq.  (1),  and  the  height  correlation  funtions 
(Ih(q)l  >,  given  by  Eq.  (2).  These  calculations  serve  to  compress  the  data  into  a 
more  useful  and  interpretable  format,  allowing  quantitative  determinations  of  the 
scaling  behavior  of  the  sputtered  surfaces.  It  is  instructive  to  first  calculate  these 
quantities  for  two  limiting  cases,  a  periodic  surface  and  a  totally  random  surface. 
The  periodic  surface,  shown  in  Fig.  9a,  is  a  high  magnification  atomic  resolution 
image  of  graphite,  collected  at  a  size  of  40  A  x  40  A.  The  random  surface 
topograph,  shown  in  Fig.  9b,  was  produced  by  substituting  the  signal  from  a 
random  noise  generator^ ^  for  the  tunneling  microscope  signal  and  allowing  the 
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STM  data  acquisition  system  to  collect  the  data  as  usual.  To  facilitate  the 
comparison,  the  random  surface  image  was  also  acquired  at  the  40  A  x  40  A  size. 
The  size  of  these  images  is  actually  irrelevant,  as  the  results  can  be  generalized  to 
images  of  any  size. 

Autocovariance  functions  for  these  limiting  surfaces  are  shown  in  Fig.  10a. 
In  general,  G  probes  the  correlation  of  surface  features  as  a  function  of  separation 
distance,  with  nonzero  values  indicating  correlated  structures.  Thus,  G(L)  takes 
on  positive  and  negative  values  and  shows  a  periodic  signature  for  the  periodic 
features  in  the  atomic  resolution  image  of  Fig.  9a.  For  the  random  surface  of  Fig. 
9b,  in  which  the  data  points  are  completely  uncorrelated,  G  *  0  for  all  lengths 
except  L  =  0.  The  large  peak  in  G  at  the  origin  is  the  variance  of  the  height  values. 
It  is  important  to  note  that  the  number  of  data  points  used  to  calculate  G(L) 
decreases  linearly  as  L  increases,  thus  there  is  greater  uncertainty  in  the  auto¬ 
covariance  function  at  large  L.  This  is  not  critical,  however,  because  the  large  L 
regime  will  not  be  used  in  the  analysis. 

The  height  correlation  function  ( I  h(q)  I  )  is  the  Fourier  transform  of  G 
and  displays  the  scaling  beharior  of  the  surfaces  as  a  function  of  w’avevector  q. 

The  sample  average  (...)  is  performed  by  rotationally  averaging  lh(q)  1 2  to  obtain 
the  radial  part  of  the  correlation  function.  As  shown  in  Fig.  10b,  the  correlation 
function  for  the  periodic  surface  clearly  show's  peaks  at  2.2  A'^  and  2.8  A-l 
superimposed  on  a  sloping  background.  These  peaks  correspond  to  the  atomic 
features  with  spadngs  of  2.24  -  2.86  A  seen  in  the  topograph  in  Fig.  9a.  As 
expected,  ( I  h(q)  1 2)  vs.  q  is  roughly  constant  for  the  random  surface  in  Fig.  9b, 
indicating  that  surface  features  are  present  at  all  wavevectors  at  approximately 
the  same  height.  Note  also  that  since  the  sampling  accuracy  of  the  Fourier 
transform  decreases  for  features  that  approach  the  image  size  in  real  space,  the 
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uncertainty  in  the  correlation  function  is  greatest  for  small  values  of  the 
wavenumber  q,  which  correspond  to  these  lengths. 

A.  Surface  Corrugation 

•  • 

Plots  of  the  corrugation  vs.  lateral  distance  for  three  2400  A  x  2400  A 
images,  Figs,  la,  2c  and  3b,  are  shown  in  Fig.  11a.  As  mentioned  previously, 
these  surfaces  make  up  a  representative  sampling  of  the  extremes  of  the  data. 

The  corrugations  for  these  surfaces  are  all  less  than  0.7  for  length  scales  dowm  to 
-  40  A,  which  indicates  structures  with  apex  angles  greater  than  110®  at  these 
length  scales.  The  topographs  in  Figs.  1-4  appear  much  rougher  than  this 
because  the  length  scale  for  the  Z  -  axis  is  highly  expanded  relative  to  the  X  -  and 
Y  -  axes.  As  a  more  careful  check  of  image  reliability  at  lengths  less  than  50  A, 
the  corrugations  for  high  magnification  images  (Figs.  5c,  6c  and  7c)  of  these  same 
representative  surfaces  are  shown  in  Fig.  11b.  In  some  cases,  this  resulted  in 
higher  values  of  o  for  these  high  magnification  images  as  the  length  scale 
approached  zero. 

Very  blunt  tips  could  also  produce  images  ^^ith  small  values  of  o  because 
the  tip  woxild  not  be  able  to  follow  sharp  features.  For  the  tips  used  to  image  rough 
surfaces  in  this  study,  however,  atomic  resolution  images  of  clean  graphite  were 
also  obtained,  indicating  the  tips  were  <  9  A  in  diameter.30  At  the  length  scales 
for  the  topographs  presented  in  Figs.  1-4,  this  tip  width  is  negligible. 

Similar  results  for  the  corrugation  were  found  for  all  surfaces  in  this  study 
except  the  high  magnification  images  of  the  surface  shown  in  Fig.  Ic.  This 
sample,  which  was  sputtered  with  the  highest  flux  (J3)  and  the  lowest  fluence 
(Q  j),  had  values  for  o  that  approached  15  for  lengths  smaller  than  10  A, 
indicating  that  these  images  most  likely  had  some  tip  artifacts  at  small  lateral 
distances. 
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B.  Auto-covariance  function 

For  the  rough  surfaces  in  this  study,  useful  information  can  be  extracted 
from  the  auto-covariance  function  G(L)  given  by  Eq.  (1)  if  certain  simplifying 
assumptions  are  made.  Because  the  mathematical  function  that  describes  the 
real  space  profile  of  a  rough  surface  is  unknown  and  difficult  to  find,  one  usually 
assumes  that  this  profile  represents  a  stochastic  process  that  can  be  described  by 
a  probability  distribution  function.24  Thus,  the  surface  profile  can  be 
characterized  by  a  knowledge  of  the  moments  of  the  distribution  function.  If  we 
assume  a  Gaussian  probability  function,  then  the  rough  surfaces  measured  here 
can  be  characterized  by  the  root-mean-square  (rms)  surface  height  and  the  1/e 
wndth  of  the  auto-covariance  function,  called  the  auto-covariance  length.  Since 
the  data  in  the  topographs  are  fitted  to  a  plane  to  remove  any  tilt  in  the  image,  the 
mean  height  of  the  surface,  <h(r))  -  0.  This  leaves  the  auto-covariance  function 
G(L)  with  beharior  for  small  L  that  can  be  approximated  as: 

G(L)  =  62  exp  (-L2//.2;.  (3) 

where  the  interface  ^^’idth  6  h  VG(0)  is  the  rms  surface  height  and  the  auto¬ 
covariance  length  /.  is  the  short  range  measure  of  the  lateral  correlation  of  the 
surface  features. 

For  the  topographs  of  sputtered  graphite,  the  auto-covariance  functions 
were  deiermined  both  parallel  and  perpendicular  to  the  fast  scanning  direction  of 
the  STM,  proriding  a  check  on  the  isotropy  of  the  images.  Significant  differences 
in  G(L)  along  these  two  directions  would  indicate  an  asymmetry  in  the  tip  or  the 
roughening  of  the  surface.  For  this  set  of  data,  however,  no  significant 
differences  were  seen  for  the  two  orientations  of  the  auto-covariance  functions. 
Therefore,  G(L)  is  shown  only  perpendicular  to  the  direction  of  scanning. 
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Table  I  gives  the  values  of  X  and  6  for  the  images  in  Figs.  1  -  3.  In  addition, 
plots  of  G(L)  vs.  L  for  some  of  the  surfaces  of  Figs.  1-3  are  shown  in  Fig.  12.  The 
initial  regions  in  these  plots  have  been  approximated  as  Gaussian  in  shape  to 
permit  the  extraction  of  an  interface  width  and  auto-covariance  length.  For  the 
lowest  fluence  Q^,  the  behavior  of  G  with  increasing  flux  is  shown  in  Figs.  12a, 
12b  and  12c,  corresponding  to  the  surfaces  in  Figs,  la,  lb  and  Ic,  respectively.  The 
fact  that  X  >  0  means  that  the  features  produced  by  the  ion  bombardment  are  not 
random,  as  one  might  expect,  but  are  correlated.  The  auto-covariance  length  X 
was  approximately  constant  for  the  two  lower  fluxes  but  decreased  to  20  A  for  the 
highest  flux,  at  which  the  G(L)  vs.  L  curve  now  looks  more  like  that  in  Fig.  10a  for 
the  random  noise  surface.  The  values  of  the  interface  width  6  for  these  curx'es 
increased  with  increasing  flux,  in  the  ratio  of  1  :  3.3  :  20,  compared  to  the  flux 
ratio  of  1  :  5  :  10.  From  Table  I,  we  obsen’e  similar  beharior  with  increasing  flux  J 
for  the  intermediate  fluence  Q2:  X  decreased  steadily  as  J  was  increased  and  the 
interface  width  6  increased,  albeit  more  slowly  than  for  the  low’est  fluence,  with 
the  ratio  1  :  2.1  :  3.4.  These  trends  in  the  parameters  X  and  6  for  the  low  and 
intermediate  Q  quantify  the  increasing  dominance  of  short  range  roughness  and 
illustrate  the  increase  in  the  random  nature  of  the  surfaces  as  the  ion  flux  is 
increased.  The  values  in  the  third  column  of  Table  I,  representing  the  highest 
fluence  Q3,  do  not  follow  the  trends  seen  for  lower  fluences.  We  see  here  that  the 
surface  with  the  largest  5  was  produced  at  the  intermediate  flux. 

Figs.  12a,  12d  and  12e  show  the  behavior  of  the  auto-covariance  function  as 
the  fluence  is  increased  at  the  lowest  flux  Jj,  and  correspond  to  the  images  in 
Figs,  la,  2a  and  3a,  respectively.  For  these  surfaces,  the  auto-covariance  length 
increased  by  -100  A  for  each  tenfold  increase  in  the  fluence  Q,  indicating  the 
formation  of  features  at  larger  separation.  The  structures  also  increased  rapidly 
in  size,  as  seen  in  the  dramatic  increase  of  interface  width  5  in  the  ratio  of  1  :  3.9  : 
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43,  for  fluences  that  increased  in  the  ratio  1 ;  10  :  100.  The  second  row  in  Table  I 
shows  similar  results  for  X  and  5  at  the  intermediate  flux  J2,  although  the 
increase  in  the  interface  width  was  more  moderate,  with  a  ratio  of  1  :  2.5  :  28.  For 
surfaces  sputtered  at  the  highest  flux  J3,  the  behavior  of  G(L)  with  increasing  Q  is 
quite  different.  The  value  of  X  was  roughly  constant  between  the  low  and 
intermediate  fluences  but  increased  abruptly  at  the  highest  Q.  The  interface 
widths  indicate  that  the  surface  is  smoothest  for  the  intermediate  value  of  Q.  This 
behavior  may  have  resulted  from  the  experimental  difficulties  in  controlling  the 
sputtering  time,  and  hence  the  fluence,  for  the  surface  etched  at  the  highest  flux 
(J3)  and  the  smallest  fluence  (Q2). 

It  is  especially  interesting  to  compare  the  rates  of  increase  of  6  discussed 
above  ^-ith  a  crude  estimate  of  the  interface  ^^^dth  which  neglects  the  presence  of 
any  smoothing  mechanisms.  In  reality,  some  smoothing  is  expected,  thus  this 
simplistic  model  should  overestimate  the  roughness  produced.  If  we  assume 
individual  ion'surface  impacts  to  be  independent  of  each  other,  we  would  expect  a 
random  surface  similar  to  that  shown  in  Fig.  9b,  for  which  X  *=  0.  As  the  fluence 
Q  is  increased,  5  would  increase  as  well,  but  the  surface  would  remain  random  in 
character.  The  arrival  of  ions,  and  hence  the  removal  of  material  from  the 
surface,  is  assumed  to  be  governed  by  Poisson  statistics,  for  which  the  standard 
deviation  in  h(r)  is  given  by  the  square  root  of  the  amount  of  material  sputtered 
from  the  surface.  The  interface  width  6  should  be  approximately  proportional  to 
this  quantity,  with  the  interplanar  spacing  as  the  proportionality  constant, 
resulting  in: 


5  .  dVQY/p, 


(4) 
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where  Y  is  the  sputtering  yield,  taken  to  be  1.5  here,  p  is  the  monolayer  area 
density  (3.8  x  10^5  atoms/cm2)  and  d  (=  3.35  A)  is  the  interlayer  spacing  for 
graphite.  The  results  of  this  estimate  for  the  three  fluences  used  in  these 
experiments  are  also  given  in  Table  I  and  have  the  ratio  1  ;  3.2  ;  10.  It  is  clear 
from  the  ratios  presented  above  that  the  roughness  increased  with  fluence  more 
rapidly  than  the  Qb'2  dependence  of  this  crude  estimate  of  6.  This  behavior  is 
quite  unexpected,  and  indicates  that  there  is  some  mechanism  in  addition  to  shot 
noise  that  is  roughening  the  bombarded  surfaces  at  high  Q.  In  other  words,  the 
individual  sputtering  events  are  not  independent  of  each  other  but  are  correlated 
in  some  way,  so  the  surface  retained  a  memory  of  previous  ion  impacts. 

Plots  of  the  auto-covariance  function  G(L)  for  the  surfaces  sputtered  at 
different  temperatures  (Fig.  4)  are  shown  in  Fig.  13.  There  was  a  marked 
decrease  in  the  interface  width  5  as  the  temperature  was  increased.  In  addition, 
the  auto-covariance  length  X  decreased  as  well,  from  100  A  to  35  A,  but  only  for 
the  highest  temperature  used  (900K).  These  values  quantify  the  decrease  of  the 
short  scale  roughness  seen  in  the  images  of  Fig.  4  but  give  little  information  about 
changes  in  long  range  correlations. 

The  auto-covariance  function  G  prorided  a  framework  for  a  quantitative 
analysis  of  rough  surfaces  in  real  space.  It  was  useful  in  demonstrating  the 
existence  of  a  surface  memory  effect,  in  extracting  values  for  the  short  range 
correlations  given  by  the  auto-covariance  length  X,  and  in  quantifjdng  the 
roughness  of  the  sputtered  surfaces  given  by  the  interface  width  6.  As  seen  in 
Table  I,  the  trends  in  these  parameters  can  be  related  to  changes  in  the  flux,  ion 
fluence  and  sample  temperature,  and  agree  with  the  qualitative  trends  observed 
in  the  topographs  of  Figs.  1  -  4.  In  particular,  the  behavior  of  the  interface  width 
as  the  fluence  Q  is  varied  indicates  the  existence  of  a  roughening  mechanism  for 
sputtered  graphite  in  addition  to  shot  noise.  The  utility  of  the  real  space  auto- 
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covariance  function  is  limited,  however,  because  tLs  issumption  of  Gaussian 
behavior  for  G(L)  vs.  L  is  followed  only  approximately  by  the  sputtered  graphite 
surfaces  and  there  is  no  detailed  theory  relating  the  observables  to  the  surface 
formation  processes. 

C.  Height  correlation  function 

A  very  useful  complementary  view  of  the  scaling  behavior  of  the  sputtered 
surfaces  was  obtained  by  using  £q.  (2)  to  calcvilate  from  the  topographs  the  height 
correlation  function  as  a  function  of  the  wavenumber  q.  The  required  sample 

9 

average  ( I  h(q)  I  )  was  obtained  by  performing  a  rotational  average  of  the  two- 
dimensional  1  h(q)  1 2  function.  For  each  value  of  q,  I  h(q)  1 2  was  sampled  at 
angular  increments  of  1®  and  averaged  over  a  range  of  180®.  Because  it  is 
symmetric,  only  half  of  the  Fourier  transform  w’as  needed  for  the  calculation. 

9 

The  variance  of  <  I  h(q)  I  )  was  also  calculated  to  provide  an  estimate  of  the 
uncertainty  in  the  correlation  function. 

Figs.  14  -  17  show  plots  of  ( I  h(q)  I  )  vs.  q  calculated  from  the  topographs  in 
Figs.  1  -  7.  The  general  behavior  of  the  correlation  functions  for  bombarded 
surfaces  is  clearly  different  from  that  shown  in  Fig.  10b  for  either  the  random  or 
the  periodic  surfaces,  again  indicating  that  the  stochastic  process  of  ion 
bombardment  has  produced  surfaces  that  are  neither  random  nor  periodic,  but 
are  correlated.  The  scaling  behavior  of  features  on  these  surfaces  can  be 
quantified  by  determining  the  correlation  length  \  and  the  functional  dependence 
of  ( I  h(q)  I  vs.  q.  The  corrugation  o  can  be.used  to  visualize  the  meaning  of  ^  in 
real  space.  As  observed  in  Fig.  11,  the  corrugation  of  each  surface  becomes  small 
for  points  separated  by  lengths  larger  than  %.  Thus,  for  such  length  scales,  the 
surface  can  be  considered  to  be  flat.  For  this  real  space  behavior,  one  would 
expect  that  the  height  correlation  function  should  be  q-independent  for  q  <  2rJ\, 
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while  it  should  decrease  with  q  for  q  >  2tU%.  This  signature  of  correlation  is 
clearly  seen  in  the  curves  of  Figs.  14  -17. 

The  creation  of  correlated  structures  by  particle  radiation  appears  to  be 
counter-intuitive.  However,  theoretical  studies  of  the  non-equilibrixim  growth  of 
interfaces20.21  indicate  that  correlation  resulting  from  random  events  is  possible. 
From  linear  response  theory  for  the  case  of  isotropic  radiation  erosion,  the 
general  form  of  the  continuum  differential  equation  describing  the  evolution  of  the 
surface  profile  as  a  function  of  time  and  position  is42,43; 

^  cc  -DgV4h  -  DyV3h.YV2h-Je(r,{h})-».Ti(r,t).  (5) 

The  first  four  terms,  which  represent  the  smoothing  mechanisms  of  surface 
diffusion,  volume  diffusion,  sputter  redeposition  and  sputter  removal  affected  by 
shadowing,  respectively,  are  illustrated  schematically  in  Fig.  18.  The  diffusion 
terms  describe  annealing  by  flow  of  material  along  the  surface  or  through  the 
bulk.  The  sputter  redeposition  term  describes  removal  of  material  by  sputtering 
and  subsequent  recondensation  at  a  different  place  on  the  surface.  The  sputter 
removal  term  incorporates  the  dependence  on  exposure  angle  6.  This  so-called 
’shadowing*  term  accounts  for  the  preferential  erosion  of  the  tops  of  prominent 
features  due  to  their  increased  exposure  to  the  incoming  flux.^^  The  final  term 
■qfr.t)  is  the  Gaussian  white  noise  for  the  incident  ions  with  a  variance 
proportional  to  the  flux  J  and  represents  the  only  surface  roughening  mechanism 
in  the  theory. 

A  simplified  linear  response  argument  was  proposed  for  radiation  erosion 
in  Ref.  44  and  compared  to  a  subset  of  our  data  in  Ref.  45.  Since  Dg  »  ,  i.e.  the 

surface  diffusion  constant  is  alw’ays  much  larger  than  that  for  the  bulk,  the 
volume  diffusion  term  was  neglected.  In  addition,  sputter  redeposition  w’as  also 
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assxuned  to  be  small  in  Ref.  45,  but  has  been  explicitly  included  here.  The 
resulting  continuum  differential  equation  is: 


=  -  DV^h  •  yv2h  -  Je(r,{h))  +  Tifr.t). 


(6) 


where  D  is  proportional  to  the  surface  diffusion  constant  and  y  contains  a  constant 
characteristic  of  the  material.  Fourier  analysis  of  Eq.  (6)  leads  to: 


3h(q,t) 

dt 


-  C3(q)  h(q,t)  +  Ti(q,t)  , 


(7) 


where  the  three  smoothing  terms  have  been  grouped  together  in  coCq),  the  healing 
rate  of  a  surface  modulation  of  wavevector  q: 

C3(q)“  Jlql +YlQ*^  +  Dlql^  .  (8) 

The  resulting  reciprocal  space  correlation  function  is 


<  1  h(q)  1 2  )  «  P  .  gjjp  2ci)(q)t)] .  (9) 

cjCq) 

According  to  Eq.  (9),  for  small  q,  ( I  h(q)  1  ^)t  is  independent  of  q,  while  for  large  q  it 
should  decrease  with  q,  indicating  two  distinct  regions  in  plots  of  (lh(q)  I  vs.  q. 
Figs.  14  -  17  show  that  this  general  behavior  is  exactly  what  was  seen  in  the 
experimentally  determined  correlation  functions. 

Fig.  14  shows  the  dependence  of  the  correlation  function  on  flux  J,  for  fixed 
fluence  Q.  The  crossover  wavevector  qo  =  2n/^  ,  which  marks  the  transition  in  the 


behavior  of  <  I  h(q)  I  from  a  constant  to  a  decreasing  finction  of  q,  is.  defined  by  .> 

y>  I  I  . 


oKqo)  1 4 1-  Eq.  (8)  predicts  thatthe  correlation  length  at  crossover,  ^  «  Q  *  Jt.  >  ' 
^assuming  %  <  and  thus  that  qo  should  remain  constant  for  the 


correlation  functions  in  Fig.  14a  and  in  Fig.  14b,  which  correspond  to  the  lower 
fluences  and  Q2.  respectively.  This  beharior  is  seen  for  the  two  lower  cur\’es  of 
Fig.  14a,  which  correspond  to  fluxes  Jj  and  J2.  In  the  upper  curve,  however, 


■which  represents  the  highest  flux  J^,  q©  has  increased.  In  Fig.  14b,  Qo  also 
increased  with  increasing  flux,  albeit  rather  slowly.  These  trends  can  also  be 
seen  from  the  values  of  calculated  from  the  crossover  wavevectors,  given  in 
Table  I.  To  investigate  the  behavior  of  the  correlation  function  for  the  limit  q  -»  0, 
the  exponential  term  in  Eq.  (9)  can  be  expanded  and  gives  the  result 
q^™Q(  I  h(q)  1 2)^  «  Jt  =  Q,  indicating  that  the  curves  in  Fig.  14a  and  Fig.  14b,  for 
which  Q  is  constant,  should  all  converge  to  the  same  point  in  the  small  q  limit. 
This  result  is  independent  of  our  choice  of  the  healing  rate  function,  o)(q),  as  long 
as  (D(q)  — ♦O  as  q  — »0.  Within  the  experimental  uncertainty,  this  is  true  for  the 
curves  corresponding  to  fluxes  Jj  and  J2  in  Fig.  14a  and  for  the  curves 
corresponding  to  fluxes  J2  and  J3  in  Fig.  14b. 

For  very  large  q,  CL)(q) «  Dq^,  and  from  Eq.  (8)  we  have  •  h(q)  I  ^)t 

J/CDq*^).  This  predicts  that  the  log/log  plots  of  ( I  h(q)  I  ^>t  shown  in  Fig.  14a  should 
be  straight  lines  with  slopes  of  •  4  for  large  enough  q  and  should  be  separated  from 
each  other  vertically  by  amounts  which  scale  linearly  -t^dth  the  fluxes,  which  are 
in  the  ratio  of  1  :  5  ;  10.  The  three  curves  all  appear  to  have  roughly  the  same 
slope,  but  the  q  dependence  is  between  q'2.5  and  q‘2-9,  and  the  vertical  separations 
of  ( I  h(q)  I  curves  at  large  q  follow  the  ratio  1:10:  400  in  Fig.  14a  and  1:10:  100 
in  Fig.  14b.  Thus,  over  the  range  of  q  in  the  plots  of  Figs.  14a  and  14b,  both  surface 
diffusion  («  q*'^)  and  sputter  redeposition  («  q*2)  contribute  significantly  to  the 
healing  function  (o(q).  The  importance  of  sputter  redeposition  has  also  been 
demonstrated  by  STM  studies  of  single  ion  impacts,46>47  which  have  shown 
considerable  transport  of  material  in  the  ■vicinity  of  impact  sites. 

The  behavior  of  the  correlation  functions  as  Q  is  varied  is  shown  in  Fig.  15. 
Although  Eq.  (9)  predicts  the  correlation  length  ^  to  be  proportional  to  the  fluence, 
we  find  from  the  data  in  Figs.  15a  and  15b,  for  the  fluxes  Jj  and  J2,  respectively, 
and  from  the  calculated  values  of  ^  in  Table  I,  that  the  correlation  length 
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increases  more  slowly  than  linearly  with  Q  for  all  fluxes  investigated.  In 
addition,  we  can  define  a  long-wavelength  interface  width  W  from  the  height 
correlation  fimctions  as  W  =  I  h(q)  I  «  (Jt)^2  =  q1/2  .  This 

asjTiiptotic  limit  of  the  linear  response  theory  represents  the  solution  of  Eq.  (6) 
with  no  smoothing  mechanisms  and  corresponds  to  the  previous  estimate  from 
Eq.  (4)  of  the  expected  interface  width  5  in  the  auto-covariance  function.  For  the 
data  in  Figs.  15a  and  15b,  Q  was  increased  in  the  ratio  1 : 10  :  100  from  the  lowest 
to  the  highest  curves.  For  a  QI'^2  dependence,  we  expect  the  interface  width  W  to 
increase  as  1  ;  3.1  :  10.  We  see  from  Table  1  that  for  the  two  lower  fluxes  and  J2, 
W  increased  in  the  ratios  1  :  9.5  :  124  and  1:11:  325,  respectively,  a  dependence 
more  like  Q  than  Q^2  ^  in  the  analysis  of  the  auto-covariance  frinction,  this 
behavior  shows  that  the  surfa.  e  roughening  from  bombardment  exceeds  even  the 
predictions  of  a  model  v.ith  only  a  stochastic  tov'’^  .ing  mechanism  included. 
This  remarkable  result,  which  \nolatcs  the  linear  response  theory  of  Ref.  44,  is 
direct  e\'idence  that  a  simple  linear  theory  cannot  explain  the  nonequilibrium 
surfaces  produced  by  the  sputtering  of  amorphous  carbon. 

For  the  highest  etching  flux,  Jg,  the  interface  width  W  increased  with  Q  in 
the  ratio  1  :  3.4  :  7.5,  which  is  closer  to  a  Ql'^2  dependence.  This  suggests  that  any 
nonlinear  roughening  mechanism  for  the  two  lower  fluxes  is  no  longer  dominant 
at  the  highest  flux.  The  energy  deposition  rate  into  the  surface  was  the  highest 
for  this  value  of  the  flux,  which  means  that  localized  heating  by  the  ion  collisions 
with  the  surface  may  be  significant.  Calculations  of  the  thermal  spike,^®  the  local 
increase  in  temperature  caused  by  a  collision  of  a  5  keV  argon  ion  with  the 
surface,  indicate  temperatures  of  -  660  K  in  the  vicinity  of  an  impact  site.  As 
already  observed,  this  temperature  is  enough  to  cause  some  enhancement  in  the 
surface  diffusion  at  the  smaller  length  scales.  Compared  to  the  lower  fluxes,  at 
the  highest  value  J3  the  impacts  with  the  surface  occur  closer  in  time  and  spatial 
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separation,  possibly  allowing  enhanced  surface  diffusion  to  effectively  smooth  out 
the  extreme  roughening  observed  at  the  lower  fluxes.  Finally,  for  large  fluences 
(ie.  long  times)  and  q  »  0,  the  exponential  term  in  Eq.  (9)  can  be  neglected, 
making  ( I  h(q)  I  ^)t  independent  of  Q.  The  experimental  results  shown  in  Figs.  15a 
and  15b,  however,  indicate  a  significant  Q  dependence  of  the  correlation  function 
at  large  wavevector  for  the  range  of  fluences  studied. 

As  a  final  test  of  Eq.  (9),  we  investigated  the  behavior  of  the  correlation 
function  { I  h(q)  I  for  sputtering  experiments  performed  at  elevated 
temperatures.  Plots  of  (I  h(q)  I  calculated  from  the  topographs  in  Fig.  4  are 
sho^m  in  Fig.  16.  These  surfaces  were  all  bombarded  at  the  intermediate  flux  and 
fluenoe  (J2  and  Q2)  but  at  progressively  higher  temperatures.  Above  600°K, 

<  1  h(q)  I  drops  more  sharply  with  q  and,  for  large  q,  has  a  tail  with  an 
approximate  q"^  dependence,  which  demonstrates  that  surface  diffusion 
dominates  the  healing  function  at  the  higher  temperatures. 

^Tiile  the  general  trends  in  the  experimental  data  are  consistent  with  Eq. 
(9),  the  quantitative  comparisons  above  rule  out  the  strict  validity  of  the  linear 
response  theory  of  Eq.  (6)  for  ion  etched  graphite  over  the  range  of  parameters  that 
we  investigated.  A  roughening  mechanism  besides  shot  noise  is  required  to 
account  for  the  Q  dependence  of  the  roughening  observed  on  sputtered  graphite. 
This  means  that  Eq.  (6)  must  have  another  positive  term,  perhaps  one  that  is 
proportional  to  h^  or  (Vh)2 

At  present,  no  non-linear  theory  exists  for  radiation  erosion,  but  a  general 
scaling  description21  has  been  developed  in  the  context  of  various  growth 
models.20  Analysis  of  the  scaling  behavior  of  the  data  may  provide  clues  to  the 
mechanism  for  surface  roughening  of  ion  bombarded  graphite.  According  to  this 
scaling  theory:^ 5 
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(Ih(q)l2)t  «  q*VF(tqZ), 


(10) 


with  F(x) «  for  small  x,  and  F(x)  -  constant  for  large  x.  This  would  predict  a 
power  law  decrease,  q*v ,  in  ( I  h(q)  1 2)^  for  q  > 

In  Figs.  14  and  15,  the  correlation  functions  calculated  from  2400  A  x  2400  A 
images  show  an  approximate  power  law  dependence  at  large  q  with  an  associated 
exponent  -v  of  order  *2.5  to  -2.9.  Because  of  the  limitations  on  the  data  density  in 
the  topographs,  however,  the  behavior  of  ( I  h(q)  I  can  only  be  plotted  up  to  q  -  0.2 
A‘l  for  these  images.  The  range  of  q  can  be  extended  to  larger  values,  however,  by 
overlapping  the  correlation  functions  calculated  for  images  of  increasingly 
smaller  scan  size.  This  was  done  in  Fig.  17  with  ( I  h(q)  I  curves  calculated 
from  the  topographs  in  Figs.  5  -  7,  which  have  sizes  of  1000  A  x  1000  A,  250  A  x  250 
A  and  50  A  X  50  A,  and  the  correlation  functions  computed  from  the 
corresponding  2400  A  x  2400  A  images.  The  plots  of  ( 1  h(q)  1  ^)t  for  these  three 
surfaces  extend  over  four  decades  in  wavenumber  q.  In  most  cases  the 
correlation  function  curves  overlap  quite  well  with  each  other.  The  discontinuity 
between  some  of  the  curves  of  differing  image  sizes  may  be  caused  by  the  fact  that 
each  image  was  leveled  before  the  power  spectrum  was  computed,  which  can 
introduce  small  vertical  scale  shifts.  This  results  in  slightly  higher  values  of 
( I  h(q)  I  )t  from  the  smaller  images.  For  the  majority  of  the  correlation  functions 
in  Fig.  17,  however,  the  curves  do  overlap  within  the  error  bars,  shown  in  Figs.  14 
and  15,  which  indicate  an  uncertainty  of  plus  or  minus  one  standard  deviation. 
More  significantly,  the  slopes  of  the  curves  at  large  q  agree  with  those  calculated 
at  smaller  q  from  the  low  magnification  images,  indicating  that  the  scaling 
behavior  observed  in  the  correlation  functions  from  large  images  (distances  >  2000 
A)  is  consistent  doT^m  to  length  scales  <  50  A,  where  the  surfaces  begin  to  appear 
smooth.  The  position  of  the  crossover  wavevector  qj,  is  approximately  the  same  in 
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curves  of  the  correlation  function  computed  from  images  of  different 
magnifications. 

Under  conditions  of  rotational  invariance,  the  exponent  z  in  Eq.  (10)  is 

related  to  v  by  z  =  2  -  a,  with  2a  =  v  -  2.  From  the  experimentally  determined 

correlation  functions,  v  -  2.5  -  2.9,  resulting  in  z  -  1.6  -  1.8  and  v/z  -  1.5.  This 

would  mean  that  for  small  q,  ( I  h(q)2)t^  -  iylz  should  increase  faster  than  linearly 

with  increasing  time  (fluence),  while  the  correlation  length  at  crossover,  ^  -  t^^z  ^ 

should  increase  more  slowly  than  linearly.  As  we  have  seen  earlier,  both 

functional  dependences  were  observed  for  the  experiments  in  this  study.  The 

values  of  a  obtained  from  numerical  studies  in  three  dimensions  are  model 

dependent;^^  for  weak  non-linearity,  a  values  in  the  range  0.15  -  0.23  were 

OA 

reported,  while  a  is  for  strong  coupling.  For  our  experiments,  we  find  a  to  be 
-  0.2  •  0.4.  The  results  embodied  in  Figs.  14  and  15  thus  appear  to  be  consistent 
with  local  growth  models,  at  least  at  the  lower  fluences.  This  demonstrates  that 
there  is  a  strong  connection  between  the  morphologies  of  surfaces  formed  by  vapor 
deposition  and  by  erosion  via  sputtering.^^ 

V.  Conclusions 

Our  STM  examinations  of  highly  oriented  p>Tol>'tic  graphite  eroded  by  ion 
beam  sputter  etching  revealed  a  number  of  dramatic  trends  in  surface  roughness 
as  the  bombardment  parameters  of  ion  flux,  fluence  and  temperature  were 
varied.  For  the  lower  fluences,  we  observed  a  steady  increase  in  roughness  as  the 
ion  flux  was  increased,  with  a  corresponding  decrease  in  the  separation 
(wavelength)  of  the  surface  features.  We  were  also  able  to  follow  the  evolution  of 
the  surface  morphology'  as  the  ion  fluence  was  increased.  Regardless  of  the  ion 
flux  used,  for  small  fluences  the  initially  flat  surfaces  showed  short  wavelength 
roughness  w-hich  evolved  into  much  larger  structures  separated  by  distances 
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which  exceeded  the  image  size  (2400  A  x  2400  A)  as  the  bombardment  was 
continued.  The  fact  that  roughness  is  produced  at  all  from  the  atomically  flat, 
cleaved  graphite  surfaces  is  significant  in  itself,  confirming  that  the  shot  noise  in 
the  incoming  ion  flux  acts  as  a  roughening  mechanism  at  the  length  scales 
investigated  here.  Because  of  enhanced  surface  diffusion,  sputtering  at  elevated 
temperatures  produced  surfaces  which  were  significantly  smoother  than  those 
bombarded  at  aunbient  conditions.  This  smoothing  effect  was  most  dramatic  for  T 
>600K 

Obser\’ations  from  the  STM  topographs  themselves  clearly  show  these 
trends,  but  only  in  a  qualitative  manner.  It  is  clearly  impossible  to  extract  much 
quantitative  information  about  the  distribution  of  length  scales  on  the  bombarded 
surfaces  by  visual  examination  of  the  topographs  alone.  In  fact,  to  the 
inexperienced  observ’er,  many  of  the  surfaces  produced  by  sputter  etching  have  a 
random  appearance.  We  therefore  relied  on  the  auto-covariance  function  G(L)  in 
real  space,  and  the  height  correlation  function  { I  h(q)  I  2)  in  reciprocal  space,  to 
pro\'ide  a  quantitative  statistical  description  of  the  scaling  and  correlation  of  the 
sputtered  surfaces.  Both  analyses  show  conclusively  that,  while  the  atomic 
ordering  of  the  cleaved  graphite  surfaces  was  obliterated  by  the  argon  ion 
bombardment,  the  resulting  surfaces  were  not  at  all  random  in  nature.  On  the 
contrary,  they  exhibited  correlated  behavnor  which  scaled  with  wavenumber  q  in  a 
well  determined  fashion.  Since  ion  bombardment  is  a  stochastic  process, 
visualized  as  random  impacts  in  time  and  space  at  these  length  scales,  there  are 
other  processes  that  are  important  in  determining  the  surface  morphology. 

The  most  significant  result  from  the  auto-covariance  function  analysis  was 
that  the  increase  in  the  interface  width  5  with  increasing  fluence  at  the  lower 
fluxes  was  greater  than  expected  from  a  simple  surface  damage  estimate, 
indicating  another  roughening  mechanism  in  addition  to  shot  noise.  At  the  lower 
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fluences,  the  decreasing  values  of  the  auto-covariance  length  X  with  increasing 
flux  suggest  that  the  surfaces  evolved  topography  with  scaling  similar  to  that  of  a 
random  surface. 

The  reciprocal  space  analysis  of  the  roughened  surfaces,  using  the  height 
correlation  function  ( I  h(q)  1 2)^,  was  the  most  useful  for  making  comparisons  to 
theories  of  surface  erosion  by  sputtering.  A  linear  response  theory  for  sputter 
erosion  was  compared  to  the  experimental  data.  This  model  incorporates  the 
smoothing  actions  of  shadowing,  sputter  redeposition  and  surface  diffusion,  and 
the  roughening  mechanism  of  shot  noise,  which  describes  the  variation  of  the 
incident  ion  flux  around  some  mean  value.  Several  of  the  predictions  of  the  model 
were  satisfied  by  the  experimental  results,  namely  that  the  correlation  function 
(I h(q) 1 2)  should  be  independent  of  q  up  to  some  crossover  wavevector  q^,  after 
which  the  curve  should  decrease  as  some  negative  power  of  q.  For  large  q  the 
theory  predicts  a  q  dependence  somewhere  between  q*2  and  q*'*.  The 
experimental  correlation  functions  for  samples  sputtered  at  ambient  temperature 
decrease  as  -q*2-'7,  which  demonstrates  that  both  sputter  redeposition  and  surface 
diffusion  contribute  to  the  surface  annealing.  For  sputtering  at  T  =  900  K,  the 
correlation  function  has  a  distinct  q*^  tail  at  large  q,  which  demostrates  that 
surface  diffusion  has  become  the  dominant  annealing  mechanism. 

However,  a  further  comparison  of  the  experimental  data  to  the  linear 
theory  revealed  several  discrepancies.  At  the  lowest  fluence  Q,  the  independence 
of  the  correlation  length  \  =  27i/qo  with  respect  to  flux  J  is  roughly  obeyed  for  the 
two  lower  fluxes,  but  not  for  the  highest  value  J3.  At  intermediate  Q,  ^  increases 
slowly  with  J.  For  constant  flux,  a  linear  dependence  of  ^  on  Q  is  expected,  but  we 
find  that  the  correlation  length  increases  more  slowly  than  linearly  with  Q  for  all 
fluxes  investigated.  The  linear  response  argument  predicts  the  long-wavelength 
interface  width  W  calculated  from  the  correlation  functions  to  be  proportional  to 
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q1/2  in  the  limit  q  — »  0.  Instead,  for  the  lower  fluxes  we  see  that  W  is  essentially 
proportional  to  Q,  indicating  the  existence  of  a  roughening  mechanism  in 
addition  to  shot  noise.  This  enhanced  roughening  behavior  suggests  that 
indiNddual  sputtering  events  are  no  longer  independent  but  are  correlated  in  some 
way.  At  higher  flux,  this  roughening  effect  is  compensated.  At  large  wavevector 
q,  linear  response  arguments  predict  ( I  h(q)  1 2)  to  be  independent  of  Q  while 
experimental  results  indicated  a  significant  fluence  dependence.  Theoretically, 
the  vertical  separation  of  the  <  I  h(q)  1 2)  curves  should  scale  linearly  with  the  flux  J 
for  fixed  Q,  but  the  separations  for  constant  Q  increase  more  quickly  than 
linearly. 

The  linear  response  theory  does  not  quantitatively  describe  the  sputtering  of 
graphite  for  the  conditions  investigated  here.  To  model  the  enhanced  roughening 
seen  at  the  lower  fluences,  Eq.  (6)  must  be  modified  with  a  new  term  proportional 
to  h^  or  {Vh)ri,  where  n  >  1.  Unfortunately,  there  is  no  analjtical  solution  for 
such  an  equation.  At  present,  we  have  no  explanation  for  the  physical 
mechanism  that  would  produce  such  a  term. 

The  fact  that  the  experimental  results  w’ere  consistent  with  a  general 
scaling  theory  applicable  to  the  evolution  of  nonequilibrium  surfaces  is 
encouraging.  The  power  law  decrease  in  the  correlation  functions  at  large  q  seen 
in  the  experimental  data  has  an  associated  exponent  of  *2.5  to  -2.9,  w’hich  is 
consistent  with  a  weak  nonlinearity.  This  behavior  extends  over  several  decades 
in  q,  as  seen  from  the  correlation  functions  calculated  from  high  magnification 
images.  In  addition,  the  scaling  theory  is  also  consistent  with  a  slower  than 
linear  increase  in  the  correlation  length  ^  with  increasing  fluence  and  a  faster 
than  linear  increase  of  the  correlation  function  with  fluence  at  small  q. 

In  general,  the  experiments  in  this  study  demonstrate  the  utility  of  the 
scanning  tunneling  microscope  as  a  tool  for  stud3ing  the  topograph}’  of  rough 
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surfaces.  The  analyses  presented  here  define  a  comprehensive  framework  for  the 
quantitative  analysis  of  surface  topography  produced  by  any  of  the  various  erosion 
and  deposition  techniques  currently  in  use.  In  this  study,  computations  of  the 
correlation  ftmctions  allowed  us  to  investigate  the  statistical  behavior  of  the 
roughness  which  developed  on  sputtered  graphite  surfaces.  It  was  then  possible 
to  make  detailed  comparisons  to  theoretical  descriptions  and  to  find  evidence  for  a 
number  of  surprising  phenomena  produced  by  sputter  erosion. 

The  exploration  of  a  larger  volume  of  the  parameter  space  of  sputtering 
variables  is  now  required.  For  example,  it  is  necessary  to  extend  the  ion  fluence  Q 
to  higher  values  to  see  if  the  enhanced  roughening  effect  continues  to  dominate. 
Also,  we  have  yet  to  determine  the  effects  on  the  sub-micron  scale  topography  of 
sputtered  graphite  of  varjnng  other  parameters  such  as  the  energy  of  the 
incoming  ions  and  the  incidence  angle.  Finally,  the  effects  of  ion  bombardment  at 
these  small  length  scales  on  materials  other  than  graphite  must  be  investigated. 
Quantitative  analysis  of  the  resulting  surface  morphologies  under  the  framework 
developed  in  this  study  will  allow  us  to  make  further  inroads  into  the 
understanding  of  the  basic  physics  of  sputter  etching  and  deposition,  and 
eventually  make  it  possible  to  predict  the  scaling  behavior  of  the  resulting 
surfaces  from  the  initial  sputtering  conditions. 
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Figure  Captious 


Fig.  1  Constant  current  STM  topographs  of  graphite  surfaces  after  sputter 
etching  with  an  ion  fluence  of  Qj  =  10^6  ions/cm^  and  fluxes  of  (a) 

6.9  X  10l3  ions/cm2  sec,  (b)  J2  =  3.5  x  ions/cin2  sec,  and  (c)  J3  = 

6.9  X  ions/cm2  sec.  The  X  and  Y  dimensions  are  2400  A. 

Fig.  2  Same  as  Fig.  1  except  the  ion  fluence  is  Q2  =  10^'^  ions/cm2. 

Fig.  3  Same  as  Fig.  1  except  the  ion  fluence  is  Q3  =  10^®  ions/cm2. 

Fig.  4  Sputtered  graphite  surface  topographs  which  show  the  dependence  of 
surface  morpholog>’  on  sample  temperature.  The  ion  fluence  and  flux 
were  fixed  at  Q2  =  lO^"^  ions/cm2  and  J2  =  3.5  x  10^^  ions/cm2  sec, 
respectively.  The  sample  temperatures  during  bombardment  were: 

(a)  T  =  300K,  Cb)  T  =  600K  and  (c)  T  *  900K. 

Fig.  5  Topographs  of  the  graphite  surface  shown  in  Fig.  la,  the  smoothest  of 
the  surfaces  sputtered  at  room  temperature,  collected  at  progressively 
higher  magnifications.  The  image  sizes  are:  (a)  1000  A  x  1000  A,  (b)  250 
A  X  250  A  and  (c)  50  A  x  50  A.  The  sample  has  been  sputtered  with 
fluence  =  lO^®  ions/cm2  and  flux  J^=  6.9  x  10^®  ions/cm2  sec. 

Fig.  6  Topographs  of  the  graphite  surface  shown  in  Fig.  2c,  a  sputtered  surface 
with  large  short  range  roughness,  collected  at  the  same  magnifications 
as  in  Fig.  5.  The  sample  has  been  sputtered  with  fluence  Q2  =  lO^^ 
ions/cm2  and  flux  J3  =  6.9  x  10^^  ions/cm2  sec. 

Fig.  7  Topographs  of  the  graphite  surface  showm  in  Fig.  3b,  a  sputtered  surface 
with  large  long  range  roughness,  collected  at  the  same  magnifications 
as  in  Fig.  5.  The  sample  has  been  sputtered  with  fluence  Q3  =  10^® 
ions/cm2  and  flux  J2  =  3.5  x  10^^  ions''cm2  sec. 
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Fig.  8 


Fig.  9 


Fig.  10 


Fig.  11 


Sawtooth  shaped  model  surface  profile  which  has  corrugation  a  =  1  at  a 
lateral  length  scale  of  L.  The  apex  angles  of  the  regular  features  are  90® 
for  these  conditions  and  no  serious  artifacts  caused  by  the  convolution  of 
tip  shape  into  the  image  should  occur  for  a  typical  probe  tip  scanning 
such  a  surface.  This  indicates  that  experimental  topographs  with  c  <  1 
should  be  free  of  tip  related  artifacts. 

(a)  High  magnification  (40  A  x  40  A)  STM  image  of  unsputtered 
graphite.  The  familiar  atomically  resolved  features  are  clearly 
seen. 

(b)  Random  surface  topograph  artificially  generated  by  feeding  the 
signal  from  a  white-noise  generator  into  the  STM  data  acquisition 
channel  and  collected  at  the  same  magnification  as  in  (a). 

(a)  Plot  of  the  autocovariance  functions  G(L)  vs.  L  for  the  surfaces 
of  Fig.  9a  (solid  line)  and  Fig.  9b  (dotted  line).  Note  the  periodic 
signature  in  G(L)  for  the  periodic  surface  and  the  lack  of  correlation 
for  all  lengths  except  L  =  0  for  the  random  surface. 

(b)  Plot  of  the  height  correlation  functions  ( I  h(q)  I  vs.  q 
calculated  from  the  topographs  shown  in  Fig.  9a  (solid  circles)  and 
Fig.  9b  (open  circles).  The  correlation  function  calculated  for  the 
periodic  surface  of  Fig.  9a  contains  peaks  which  correspond  to  the 
atomically  resolved  features  seen  in  the  STM  topograph.  On  the 
other  hand,  the  approximately  constant  value  of  ( I  h(q)  I  for  the 
random  surface  indicates  features  of  roughly  equal  height  at  all 
wavevectors. 

Plots  of  the  corrugation  a  vs.  lateral  length  calculated  from: 

(a)  the  2400  A  x  2400  A  topographs  of  sputtered  graphite  shown  in  Fig. 
la  Garge  dashes),  Fig.  2c  (solid  line)  and  Fig.  3b  (small  dashes). 
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(b)  the  high  magnification  (50  A  x  50  A)  topographs  of  the  same  three 
surfaces  as  in  (a),  which  are  shown  in  Fig.  5c  Oarge  dashes),  Fig. 

6c  (solid  line)  and  Fig.  7c  (small  dashes). 

These  surfaces  represent  the  extreme  limits  of  roughness  in  the  data 
set.  Note  that  o  <  1  for  all  length  scales,  indicating  the  absence  of  tip 
related  artifacts. 

Fig.  12  Plots  of  autocovariance  functions  that  show  the  changes  in  the  behavior 
of  G(L)  vs.  L  for  increasing  flux  J  (a,  b  and  c:  calculated  from 
topographs  in  Figs,  la,  lb  and  Ic,  respectively)  and  for  increasing 
fluence  Q  (a,  d  and  e:  calculated  from  topographs  in  Figs,  la,  2a  and  3a, 
respectively).  The  initial  regions  have  been  approximated  as  Gaussian 
in  shape  to  permit  the  extraction  of  the  auto-covariance  length  X  and 
the  interface  width  5. 

Fig.  13  Temperature  dependence  of  the  autocovariance  function,  with  ion 

fluence  and  flux  fixed  at  Q2  =  101*7  ions/cm^  and  J2  *  3.5  x  lO^^  ion&'cm^ 
sec,  respectively.  G(L)  is  calculated  from  the  topograph  in: 

(a)  Fig.  4a,  T  =  300K,  (b)  Fig.  4b,  T  =  600K,  (c)  Fig.  4c,  T  =  900K 
Note  the  decrease  in  the  interface  width  6  =  VG(0)  .  The 
autocovariance  length  A  decreases  significantly  only  at  the  highest 
temperature,  T  =  900K. 

Fig.  14  Plots  of  height  correlation  function  vs.  wavevector  q  showing  the 

dependence  on  flux  J  of  the  correlation  function  ( I  h(q)  1 2),  the  interface 
width  W  and  the  correlation  length  %  for  the  two  fluences  (a)  Q ^  =  IOI6 
ions/cm2  and  (b)  Q2  =  10 1^  ions/cm^.  The  curves  corresponding  to  the 
three  fluxes  are:  6.9  x  10l3  ions/cm^  sec  (filled  circles),  J2  =  3.5  x 

IOI4  ions/cm2  sec  (open  squares)  and  Jg  =  6.9  x  lOl'^  ions/cm^  sec  (filled 
triangles).  The  values  of  the  crossover  wavevector,  Qo  =  2tJ\,  are 
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indicated  in  each  curve  and  a  l/q2-7  dependence  for  large  q  is  shown  for 
comparison  in  (a).  The  error  bars  represent  plus  and  minus  one 
standard  deviation. 

Fig.  15  Dependence  of  the  correlation  function  ( I  h(q)  1 2),  the  interface  width  W 
and  the  correlation  length  ^  on  fluence  Q  for  the  two  fluxes  (a)  Jj*  6.9  x 
10^3  ons/cm2  sec  and  (b)  J2  =  3.5  x  10^^  ions/cm^  sec.  The  curves 
corresponding  to  the  three  fluences  are:  Qj  »  10^6  ions/cm2  (open 
circles),  Q2  *  ions/cm^  (filled  squares)  and  Qg  =  10^®  ions/cm2 
(open  diamonds).  The  error  bars  have  the  same  meaning  as  those  in 

Fig.  16  Dependence  of  the  correlation  function  ( I  h(q)  1 2)  on  sample  temperature 
T  during  bombardment:  T  =  300K  (open  squares),  T  =  600K  (filled  circles) 
and  T  =  900K  (filled  diamonds).  A  1/q^  dependence  is  shown  at  large  q 
for  comparison.  The  smoothing  of  the  surfaces  sputtered  at  higher 
temperatures  is  indicated  by  the  overall  decrease  in  the  values  of 
( I  h(q)  1 2),  Note  the  q”^  tail  at  large  q  for  sputtering  at  T  =  900K. 

Fig.  17  Plots  of  height  correlation  functions  shoeing  the  consistency  of  the 
scaling  beharior  of  ( I  h(q)  1 2)  over  approximately  four  decades  in 
wavevector  q.  The  correlation  functions  were  calculated  from 
topographs  of  sizes  2400  A  x  2400  A  (filled  triangles),  1000  A  x  1000  A 
(open  squares),  250  A  x  250  A  (filled  diamonds)  and  50  A  x  50  A  (open 
circles). 

(a)  <  I  h(q)  1 2)  calculated  from  topographs  of  Figs,  la,  5a,  5b  and  5c  for 
sputtering  conditions  =  10^6  ions/cm2  and  6.9  x  10^3 
ions/cm2  sec. 

(b)  <  I  h(q)  1 2)  calculated  from  topographs  of  Figs.  2c,  6a,  6b  and  6c  for 
sputtering  conditions  Q2  =  10^7  ions/cm2  and  Jg  =  6.9  x  10^^ 
ions/cm2  sec. 
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(c)  <  I  h(q)  1 2)  calculated  from  topographs  of  Figs.  3b,  7a,  7b  and  7c 
for  sputtering  conditions  Q3  =  10^8  ions/cm2  and  J2  -  3.5  x  10^^ 
ions/cm2  sec. 

The  curves  do  not  overlap  exactly  because  of  small  vertical  shifts  in 
{ I  h(q)  1 2)  produced  by  the  leveling  of  the  different  size(^topographs._  .Note 
the  consistency  of  the  slopes  for  each  set  of  curves.  \ 

Fig.  18  Schematic  diagram  of  the  four  smoothing  mechanism^^s^  the^* 

corresponding  terms  in  Eqn.  (5),  the  general  continuum  differential 
equation  for  the  surface  profile  derived  from  linear  response  theory. 


Table  I.  Summary  of  characteristic  lengths  (in  A)  calculated  from  the  STM 

topographs  of  Tigs.  1  •  3.  The  autocovariance  length  X  and  the  interface 
undth  5  were  determined  from  plots  of  the  autocovariance  function  G(L) 
vs.  L.  The  correlation  length  \  and  the  interface  width  W  were 
determined  from  plots  of  the  height  correlation  function  (Ih(q)l  2) 
vs.  q.  In  addition,  the  bottom  row  shows  values  for  the  interface  width 
5poicson  calculated  from  an  estimate  of  the  surface  roughening  due 
solely  to  shot  noise.  For  this  estimate,  the  roughening  follows  the 
Poisson  statistics  of  the  random  ion  bombardment  (and  material 
removal)  process. 
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Summary  of  Characteristic  Lengths  (in  A)  Calculated 
from  STM  Topographs 
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